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Abstract 

(N ; 

^ . We calculate the ultra-relativistic Bose-Einstein condensation temperature of a 

complex scalar field with weak A(<I>^<^)^ interaction. We show that at high temper- 
CSI ' ature and finite density we can use dimensional reduction to produce an effective 

. three-dimensional theory which then requires non-perturbative analysis. For sim- 

I plicity and ease of implementation we illustrate this process with the linear delta 

' expansion. 

p_; 1 Introduction 

(D : 

I The inclusion of finite densities of conserved charges in thermal field theory poses well 
I known problems in the study of phase transitions. On one hand, perturbative analyses 
^ ■ of charged scalar fields [1, 2, 3] give information about the phase structure, but cannot 
■ probe the critical point. Standard perturbative calculations are plagued with infrared 
divergences and after dealing with these, the asymptotic expansion breaks down. On the 
other hand, finite charges cannot be easily represented in lattice Monte-Carlo simulations. 
Encoding a non-zero charge in the grand canonical ensemble renders the action complex 
and therefore useless as a statistical weighting [4]. These obstacles invoke a demand for 
alternative methods. 

Several authors have recently considered the effect of repulsive interactions on the 
condensation temperature of a dilute Bose gas, a non-relativistic problem. It was pointed 
out by [5] that the leading correction could be isolated in the static Matsubara mode and 
a mean field calculation could be performed on this mode alone. Since then the static 
mode has been studied using the expansion [6, 7], with the linear delta expansion 
[8], and now numerically [9, 10, 11]. 

In this paper we consider the effect of interactions on the transition temperature for 
scalars in the ultra-relativistic limit. By this we mean that we are at the high temperature 
limit, but that the typical self-energy corrections, S, the chemical potential, /z, or the 
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cube root of the charge density, p, can be of the same order T E ~ ~ p^/^. This 
hmit is appropriate for the study of high temperature symmetry breaking where Bose- 
Einstein condensation and spontaneous symmetry breaking have interesting similarities. 
The relevant finite temperature and density 4D Feynman diagrams are difficult to handle 
and only Jones and Parkin [12] include the setting-sun diagram for the self-energy. To 
avoid these problems, we use dimensional reduction [13, 14, 15, 16] and this is the main 
focus of our paper. In our regime, we will show that dimensional reduction for a relativistic 
model is much more complicated than in the non-relativistic limit of our model which 
studied in [5, 6, 7, 8, 9, 10, 11], but we will show it is still manageable. As in the 
non-rclativistic case, we take advantage of the good infrared behaviour of the non-static 
Matsubara modes. This allows us to integrate these modes perturbatively in an attempt 
to simplify our calculation. 

Following dimensional reduction, the problem of relativistic Bose-Einstein condensa- 
tion has been reduced to the study of a phase transition in an effective three-dimensional 
theory at zero temperature and zero density. Standard non-perturbative methods can 
then be employed to study this model but we do not investigate the many alternatives 
here, merely choosing one for exemplary purposes. Here we use the linear delta expansion 
for its simplicity and resemblance to standard perturbation theory. 

By way of contrast, the only other study of a relativistic finite density system using 
dimensional reduction we know of is a study of QCD [17], rather than the Higgs sector 
studied here. The resulting 3D action is complex, making subsequent numerical analysis 
of [17] more complicated. 



2 Dimensional Reduction 

We begin by considering a relativistic system of bosons described by a complex scalar 
field theory. We encode a conserved charge by working in the grand canonical ensemble 
with chemical potential p. The partition function is then 

Z = j[d^^][d^]&q){-S} (1) 
where field integrations are periodic over imaginary time (3 = 1/T and the action is given 

by 

The charge density is obtained from Z as follows: 

T 1 dz 

where we let the volume V tend to infinity. 

The periodicity of the fields is made explicit by a mode expansion and the non- 
static modes are integrated perturbatively^. This not only gives an overall factor to Z 
but also renormahzes the parameters of the static mode, The result is an effective 
three-dimensional theory whose dependence on the temperature and chemical potential 
is contained within its mass and coupling: 

-'^Alternatively, one can calculate static quantities in the 3D effective theory and in the full 4D theory. 
By matching these results one can then relate the coefficients of the 3D effective theory to those of the 
4D one [16]. However there is a lack of results in the literature for our model at high densities. 
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Figure 1: Contributions to the mass of the static mode. Bold hues correspond to non- 
static internal lines, thin lines correspond to the static mode. The last diagram contains 
the one-loop vertex counterterm. 



Z = exp{(3VF{ii)} J exp{-53D}. (4) 

Writing F — F^^) + XF^^^ + 0{\^), the quadratic part of the action can be integrated 
to give 

F^'\pi)^Tj: [-^^lnA{u;^,p). (5) 
and, since we are dealing with bosons, a;„ ~ 



where A(a;„,p) = [{un — i^J^)^ + + 
27mT. The leading perturbative correction is given by the figure-of-eight diagram: 



(6) 



We use dimensional regularization in the MS scheme, making the replacement 

d^p r (e<M^Y f 



(27r) 



An 



{2n: 



p 

3-2e 



(7) 



and subtracting only the terms which are divergent as e — > 0. M is an arbitrary renor- 
malization scale and 7 is the Euler-Mascheroni constant. 

The ultra-relativistic limit is defined by p ^ or, equivalently, T ^ m. In order 
to avoid confusion over the different expansion parameters we will set m = 0. Since 
we shall handle the infrared region non-perturbatively, this does not create additional 
problems. We keep corrections of 0(//^/T^) as we will find these are in fact 0{\) for the 
critical theory. Ignoring any //-independent terms, the factor F is then given by (see [2] 
for details) 
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+ 



A 



47r2T2 47r2 



+ O (A^) 



(8) 



where by C>(A2) we mean 0{ii'^/T\ Xfi^T^, A^). 

The non-static modes also give corrections to the static mode action. We write the 
effective 3D action as 



= J d^x [V^'^V* + r$t$ + , 



(9) 



where we neglect the higher dimensional operators ($^$)" for n > 3. This is firstly on the 
grounds that their coefficients are high order in A and further that all non-renormalizable 
interactions are suppressed by factors of (for further discussion of this point see 

[16]). It can be checked that in the calculation which follows, these neglected terms 
contribute at higher order in the A-expansion than is considered. 



The fields should be scaled by vT in order to remove an overall factor of 1/T from the 
action and to the order we work they undergo no field renormalization from the non-static 
mode integration. 

Referring to figure 1, the 3D mass is given by 



r = -/x' + Si + S2 + S3 + E4 + S5 + 0{\^ 



where 



El = / A(a;,,p) 



27r2T2 



+ o (a^) 



A^T^ 



1 Me^ 

- + 21n^ + l-,+ ^ + 0(A) 



(10) 

(11) 

(12) 
(13) 

(14) 



Though S4 involves a static internal line it should be included since it will not arise 
from the effective static 3D theory^. This contribution is in fact well behaved in the 
infrared and to leading order we find (see appendix) 



S3 + E4 = o + e)(A''r^). 



(15) 



Finally we include the diagram with the one-loop vertex counterterm which is given 
by (5A2r/87r2e)($t$)2. This is 



, SA^T 1 „ /• ^ ^ 



5A^ 
127r2 



1 1 Me^ ^ C/(-l) 



(16) 
(17) 



We may now sum all the contributions to equation (10) giving 

AT2 3/^2 (I , l,„Me^ , ^ _ , C/(-l)~ 



1 1, .^/x2\ 



27r2T2 47r2 \^2e ' 3 47rT ' * ' ' C(- 
We shall find the coupling to be given with sufficient accuracy by 

M = AT + 0(A2). 

The coupling is now dimensionful due to the scaling of the static fields by \/T. 



(18) 



(19) 



^Notc that at this order, integrating out the heavy modes actually leaves one with a non-local action 
with terms such as / cfixtfix' 

|$(x)pB(x,a;')|$(a;')P where B is an O(A^) bubble diagram. It is when 
approximating the 3D theory by the purely local one (9) that contributions, such as S4, coming purely 
from non-local terms in the exact non-local effective theory must not be forgotten. For instance, the 
same type of non-local B term should also lead to a diagram similar to S2 but with a light petal on 
the top and one has to check that this is of lower order than required. Such problems suggest that the 
matching of Green functions approach to dimensional reduction [16] might be simpler. 
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Use of equations (3) and (4) with (8), (9), (18) and (19) gives 



OF 
dfj, 



2fj,T 



X 

2 A 



1 + 
/i 



OX 



+ 



27r2T2 47r2 



+ O (A^) 



(20) 

($t$), (21) 



where ($'''$) denotes the Green function evaluated in the effective 3D theory. This we 
cannot calculate perturbatively since the expansion will break down when probing large 
length scales greater than 1/u. 



3 Linear Delta Expansion 

To evaluate the quantity ($^$) we need a non-perturbative method as the 3D sector 
retains all the infrared divergences of the 4D theory. The effective 3D theory is studied 
at zero temperature and density so the problem is greatly simplified. Any standard 
non-perturbative method can be used at this point. 

We will use LDE, the Linear Delta Expansion, though the method is also known by 
several other names (see [20] for a brief summary). LDE has been used successfully in 
many situations, including studies of scalar theories at non-zero density such as [8, 18, 
19, 12, 20]. In toy models, where exact results are achievable, LDE is known to produce 
convergent results and to do so much faster than alternatives, for instance see [21, 22] 
and references therein. In full QFT, LDE has also often proved to be better than other 
methods [23]. 

We begin by defining Ss which interpolates between our 3D action (from which we 
now drop the subscript) and some soluble action Sq as 5 varies from 1 to 0: 



S ^ Ss = SS+{l- 6)So. 



We are free to choose 
such that 



Ss = J (fx [V$^V$ + iSf - Sn^ + 5r)$^$ + 5m($^$)^ 



(22) 
(23) 
(24) 



Any physical quantity V is evaluated as a power series in S to some finite order. 
This quantity will generally have some dependence on Q which we fix by some specified 
criterion which we take to be the principal of minimal sensitivity (PMS) [24] : 



= 0. 



(25) 



<5=i,n=n 



This variational procedure allows for the emergence of non-perturbative behaviour. 
The Green function we require can now be written as 



d^p 
(2^ 



where 



Gs{p) = P^ + n^- Sn^ + Sr + Es{p) 



(26) 
(27) 
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Figure 2: Contributions to the self-energy for the effective 3D theory. The filled dot 
corresponds to a 5(r — Q^) mass insertion. 



The criterion defining the transition temperature is that the correlation length in the 



original theory be infinite. This can be expressed as ^(0) 
imposing 

Sr + E5(0) = 0. 
Use of this relation in equation (26) gives 



5=1 



which is satisfied by 
(28) 



and expanding to second order in S we have 
(Pp 1 



(27r)3 p2 + 



1 + o . + 



S,(p)-S,(0) 

p2 + n2 



(29) 



+ 0{5^). (30) 



The only momentum dependent contribution to the self energy results from the setting 
sun diagram (see figure 2): 



d^k d^q 



1 



1 



1 



(27r)3 (27r)3 (A;2 + ^2) (^2 + ^2) ((A; + g + p)2 + J^2) • 



(31) 



We continue our use of dimensional regularization in the MS scheme, using the same scale 
M as in our 4D heavy mode calculations. The required integrals (see, for example, [8] 
and [25]) are 



L 



p p2 + 02 

Ef(0) 

^3D 



An 
86^u^ r 1 



{Any 



4e 



In 



M 
M 1 



l + 2.(ln- + l)+C)(.^) 
+ 0{e 



3Q 



p p2 + Q2 



(p) = 



(4^[2; + ^^^2Q + l-2^^2 + ^(^) 



(32) 
(33) 
(34) 



Further integrals may be derived from (32) by successive differentiations with respect to 
Summing the contributions to the Green function, we find 



45^1^2 4 
In -. 



An 2(47r) 8(47r) (47r)3Q 3 



(35) 



The divergences cancel and there is no need to invoke any counter terms. We now apply 
the PMS condition to ($^$) and arrive at 



^ u /2 A 



. , In 
TT V3 3 
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(36) 



The question of which solution to choose can be answered by comparison with the solution 
in the large- A?" case for — > 2. Referring to [6] and [7] we find that the positive solution 
is appropriate, giving 



(47r)2' 

61n- 
3. 



1.314. 



(37) 
(38) 



Inserting this into equation (21) we have for the critical density 



2^+(4^(^-^^>^M 



(39) 



The chemical potential is an unwanted free variable in this expression and we must 
use relation (28) to constrain at the transition temperature. Along with the setting 
sun contribution to the self energy, we also have to order S'^ (see figure 2) 



3D 



46u 



pp2 + n2' 



^3D 



45' 



''u{^f -r) [ 
Jp 



1 



(p2 + n2)2- 

Evaluating and summing these contributions gives 



(40) 
(41) 
(42) 



= 0) 



5un 8(52m2 



TT (47r)2 
which upon substituting into (28) gives 



1 , M 1 
4^ + ^"3Q-2. 



+ 



5'^ur 



2n 



27rQ 



+ 0{6^), (43) 



SuD, S'^ufl S'^u 
Sr = z h — 



2,, 2 r 



TT 



2tt 



27r2 



- , M 



(44) 



We apply the PMS condition to r, giving ft — u/tt, and insert this value into equation 
(44): 



u 

2^ 



1 , M 

— + In \- Cr 

4e u 



, TT 3 

In - + - 
3 2 



1.546. 



(45) 
(46) 



The dependence of r on the scale M is exact because of superrenormalisability and 
agrees with that noted elsewhere (e.g. [11]) even though it appears here in the context of 
a particular non-perturbative calculation. Comparing with equation (18) we have 



1 - 



3/i 



2 3A 
+ 



1 



27r2r2 47r2 V 2e 



1 , Me^ 

- In 

3 AttT 



-7 + ^)+0(A^) 
A^T^ r 1 , M 



(47) 
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The divergences cancel and we can obtain in terms of the critical temperature: 



^^\T 



1 + 



87r2 



Or 



T \ 
In — + 77mu In A + a^^ j + 0{\^) 



6, 



C(i,-i) 
C(-i) 



- 6c^ PS 3.270 - 6c^ 



(48) 
(49) 

-6.007. (50) 



It is now clear that for the critical theory, //^/T^ ~ ^i^- Substituting into the equation 
for the charge density gives 



3v^ 

Oixho = 5, 



1 + 



A 

8^ 



Vrho 
Omu + 77-3/ 



"rfeoln ^ + 77rfeoln(A) + a^ho^ + C(A^ 
6, 

3.9236 - 6cr - 3/ -21.759. 



(51) 
(52) 
(53) 



Equation (51) relates the critical density to the critical temperature for the Bose-Einstein 
condensation of a complex scalar field in the ultra-relativistic limit. Though (51) looks 
like an expansion in A we stress that the result is non-perturbative because of the severe 
IR problems in calculating ($^$) and self-energies at the critical point in the three- 
dimensional theory. We should expect p to have some renormalization scale dependence 
since we have included a one-loop vertex counterterm in the calculation. 



4 Discussion 

The leading term in our expressions for the critical chemical potential (48) or cquivalently 
for the critical density (51), are the usual leading high temperature results, /iq = ^X/3T 
and po = ^JX/27T^. 

The A^ In (T/M) term is exactly that expected from the running of the 4D coupling A 
in the leading term using perturbation theory where one finds 

A(M.)=A(MO + |^ln(^^j. (54) 

The perturbative result is appropriate as the leading behaviour comes only from the 
heavy modes, and these are dealt with perturbatively in this calculation. Thus we find 
that once the implicit scale dependence of A is accounted for, our results for p and p are 
actually independent of the scale M as all exact physics results should be. 

The A^ln(A) term comes directly from our expression (45) for r and in particular 
comes from dependence on the scale M which is exact for for the super renormalisable 
theory. In the context of a Green function matching approach, as discussed for instance 
in [16], this term might be described as running the 3D mass r from the scale M ~ T, 
used when dealing with the heavy modes, down to w = AT appropriate for the 3D static 
theory. Overall then our expressions for the critical chemical potential and density agree 
with our expectations from other calculations. 

In terms of actual numbers, for A = 1/8 the fractional corrections to this leading 
term are of the order of a few percent, e.g. in units of the scale M, for T = 10.0, the 
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-0.015 ; 

-0.02 ; 
-0.025 ; 
-0.03 ; 

-0.035 



-0.045 ; 

Figure 3: Fractional fi correction, fi/fio " 1, against the critical T in units of M for 
A = 1/8. Dots are data points from the full 4D Jones and Parkin method with best fit 
curve connecting them. Other line is the dimensional reduction form with LDE used for 
3D non-perturbative calculation. 

critical chemical potential is fx — Ho x (1 — 0.010), /Xq = 2.04 and the critical density is 
p = po X (1 - 0.058), po = 68.0. 

A good test of our central result, the dimensional reduction, is to compare our formulae 
for the critical chemical potential with that extracted from the results of Jones and Parkin 
[12]. They also use the linear delta expansion but they apply it directly to the full four- 
dimensional theory anywhere in the symmetric phase. They can in principle study any 
temperature and density, though they limit their analysis to high temperatures. We use 
dimensional reduction so our method is always limited to high temperatures. However, 
this brings several benefits to us as further fields, including fermions, can be added with 
great ease and our non-perturbative effort is much less. The effective three-dimensional 
theory derived here can in principle be studied using any non-perturbative technique, 
including Monte-Carlo since the effective action is real. While we have only looked at 
the critical point, we have also given the critical density which is the actual measured 
quantity. Finally, our results are completely analytic, while Jones and Parkin can only 
find numerical solutions to their equations, though the numerics are relatively straight 
forward. 

Turning to the details of the results, we find that the Jones-Parkin method gives the 
same qualitative behaviour as our results for A = 1/8, T/M = 0.5 . . . 10.0. The fractional 
correction to the leading p behaviour, Ap := fi/ fiQ — l is shown in figure 3. However their 
results for these parameters [28] are best fitted with slightly different coefficients^, namely 
oirho,jp/oirho = 1-2(1), and choosing ■qrho,jp = Vrho we find ar/io,jp-F77r/io,jpln A = -17.8(2) 
compared with our value —18.5. In our calculation the a-aiu value was set by perturbation 
theory and can not be exact. Since Jones and Parkin use an entirely non-perturbative 
method, even though both methods ought to be valid in the parameter range considered, 
there are hkely to be small differences between the two results. These ought to be 0(A) 
fractional corrections so the results seem to be consistent as far as they go. 

Finally, we can illustrate the power of our work by noting that any non-perturbative 
calculation can be used with our formulae*^ for dimensional reduction in the presence of 
a large density. For instance, since this work was completed, two lattice Monte Carlo 

•^A shift from MS to MS scales is needed. 

^Note that the normalisation of fields and the definition of u or A often differ from those used here 
by simple constant factors. 
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studies of the 3D action relevant here have appeared [9, 10, 11]. While these numerical 
results were produced for the non-relativistic 4D Bose-Einstein condensation problem, 
we can just as easily apply them to our relativistic case. The two studies are completely 
consistent but for definiteness we use the values in [11] which gives c^^iat = 0.671(1) and 
flat = 0.57(1). In this case, the a and r] terms of both and p expressions are fixed by 
dimensional reduction so the only difference is in the a parameters, and we find that the 
lattice data leads to ajnu,iat — —0.756 and a,.ho,iat ~ —13.2 compared to our values of —6.01 
and —21.8 respectively. While these appear to be large differences, in physical quantities 
such as the critical density p the constant coming from a is in fact overwhelmed by the 
contribution from the In A term for A <^ 1 which is where the dimensional reduction 
process is vahd. The difference is more of a comment on the efficacy of different non- 
perturbative approximations in 3D calculations (see [10] for a good comparison) than 
particularly important to our results. 

In conclusion, wc have shown that by organising the modes into those which can be 
pertubatively integrated and those which cannot, we minimise the non-perturbative effort 
needed to study Bose-Einstein condensation at relativistic temperatures and densities. 
The method is economical and provides a reliable estimate of the critical temperature 
and critical density. 
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Appendix 

In this appendix wc consider how the sunset diagram can be decomposed into contribu- 
tions from different modes. 

We begin by considering the /i = case and state the results (see, for example, 
[26, 27]) 



E 



1 1 1 



Jpi,P2 (^ni +Pl + m^) (^na + Pl + "^^) ((^m + ^na)^ + {Pl + V2Y + m^) 



711,712 = — 00 t'^'t''' V^ni ' /'i ' / ^^712 

1 

and 



1 , M 1 

h In h - 

4e 3m 2 



O [ (A.1) 



Pi,P2 (Pi + m2) (p2 + rri^) ((p^ + p^)2 + ^2) (^^^y 



1 M 1 



(A.2) 



We shall split (A.l) into purely non-static internal lines {rii 7^ 0;n2 7^ 0;ni 7^ n^)-, one 
static internal line, (ni 7^ 0; ^2 7^ 0; rii — n2 and ni 7^ 0; ^2 = and rii = 0; ^2 7^ 0), and 
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purely static internal lines, (ni 
write 



0;n2 = 0). This covers all the possibihties and we may 



ni,n2=— oo 



m 7^712 



ni E 

n2>0 



(A.3) 



Denoting these contributions to (A.l) as /nonstatio -^mixed and /static; we immediately see 
that /static is given by equation (A. 2). 

The contribution with one static line can be calculated by taking the static line to 
be massless. This does not cause infrared divergences and is appropriate for the critical 
theory. We may also set m = in the other propagators since T 3> m. We thus have 



' mixed 



'Pi,P2 pi pI + {pi + P2Y + 

'e^W^y r-^(l/2 + e) 




1 



, ITT- 

+ o{- 



1 



(47r)3/2-cr(l + 2e) Jp, pi (pj + ujiy/^~' ^ [f 



r2(l/2 + e) 



2e(e- l/2)(47r)3-2^ 



+ 



3C(46)r2(l/2 + 6) 
e(e-l/2)(47r)3-2^ 



+ - 



1 , M 1 

h In h - 

4e 2T 2 



(A.4) 

(A.5) 

(A.6) 

(A.7) 
(A.8) 



Finally, the purely non-static contribution can be found by subtracting the other 
contributions from (A.l): 



' nonstatic 



{Any 



1 , M l' 

h In h - 

.4e 2T 2. 



(A.9) 



Turning to the case where 7^ 0, /nonstatic is unchanged at leading order in the 
high temperature expansion since ^ <S T and we may clearly take the — > limit 
without causing infrared divergences, /mixed is also unchanged upon choosing the static 
line to be critical. Up to corrections of 0{fi'^/T'^) and 0{m/T), the non-static and mixed 
contributions to the sunset diagram sum to zero. 
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